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ABSTRACT 

In order to determine the column density of a component of an ice from 
its infrared absorption features, the strengths of these features must be known. 

The peak positions, widths, profiles, and strengths of a certain ice component’s 
infrared absorption features are affected be the overall composition of the ice. 

Many satellites within the solar system have surfaces that are dominated by H 2 O 
or N 2 and ices in the interstellar medium (ISM) are primarily composed of H 2 0. 

The experiments presented here focus on the near -infrared absorption features 
of CO, C0 2 , CH 4 , and NH 3 (^=10,000-4,000 cm -1 , A=l-2.5/ixm) and the effects 
of diluting these molecules in N 2 or H 2 0 ice (mixture ratio of 5:1). This is a 
continuation of previous results published by our research group. 

Subject headings: Ices: IR spectroscopy - Satellites: general - Satellites: surfaces 
- Spectroscopy 


1. Introduction 

A variety of frozen volatile species ( “ices” ) are observed in the mantles coating the dust 
grains of the dense interstellar medium (ISM), where they are dominated by H 2 0 or N 2 (e.g., 
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Gibb et al. 2004). and on the surfaces of planetary bodies in the outer Solar System (e.g.. 
Roush 2001). Since the spectral properties (peak positions, spectral profile structures, and 
absorption strengths) of most species may be significantly altered by dilution, it is important 
to study them in an appropriate ice matrix. Here, we extend our previous study of the near- 
infrared spectral properties of pure ices (Gerakines et al. 2005) to those of molecules present 
within H 2 0- and N 2 -dominated ice mixtures. 

The compositions of interstellar icy grain mantles are dominated by H 2 0, by about a 
factor of 2 over the next-most abundant ice species, with minor components dependent upon 
the specific environment probed along a chosen line of sight. Some of the most significant 
minor components of interstellar ices include CO (1-10% relative to H 2 0), C0 2 (10-30%), 
CH4 (1-5%), and NH 3 (~1%; with some dependence on the environment sampled in the 
observed line of sight (e.g., Gibb et al. 2004). Moreover, infrared spectroscopy of the vi- 
brational absorption features of these minor components indicates that while H 2 O-poor ice 
components do exist, a significant majority of icy mantle species appear to be diluted within 
an H 2 0-ice matrix (e.g., Tielens et al. 1991; Chiar et al. 1995; Boogert et al. 1996; Gerakines 
et al. 1999; Oliveira et al. 2009). 

In the outer solar system, many cold planetary surfaces appear or are thought to be 
dominated by II 2 0 (see, e.g., Hudson & Moore 2001), such as those of cometary nuclei 
(Roush 2001; Gibb et al. 2003), some Galilean satellites (e.g., McCord et al. 1998; Hibbitts 
et al. 2000; Hansen & McCord 2008), many of the Saturnian satellites (e.g., Filacchione et al. 
2007; Lopes et al. 2008), and many Kuiper Belt Objects (Jewitt & Luu 2004; Trujillo et al. 
2005). Spectral modeling of the reflectance spectra indicate trace species on these surfaces 
such as CH4, C0 2 , and NH 3 (e.g., Hibbitts et al. 2003; Cruikshank 2007; Schaller & Brown 
2007; Clark et al. 2008). Molecules diluted in N 2 have been identified on planetary bodies 
such as Pluto and Triton (e.g., Cruikshank et al. 1998; Quirico et al. 1999; Buie &; Grundy 
2000; Grundy et al. 2003) and N 2 and CH 4 have been detected on large Kuiper Belt Objects 
(KBOs) such as Sedna (K130) (Barucci et al. 2005) and Eris (Brown 2008). 

By and large, the current state of knowledge about the composition of interstellar ices 
comes from observations and laboratory simulations of spectra obtained at mid-IR wave- 
lengths (A = 2.5 — 25 pm). This is due to the fact that the strongest absorption features of 
most small molecules in the solid phase fall in the mid-IR spectral region. These features 
are usually caused by fundamental molecular vibrational modes and are therefore the most 
sensitive tracers of ice composition (e.g., the 3.05 pm O-H stretching vibration feature of 
H 2 0 or the 4.26 pm C=0 stretching vibration feature of C0 2 ). 

In contrast, the majority of what is known about Solar System ices has been obtained 
from the spectra of reflected sunlight in the near-IR spectral region, where the absorption 
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features due to overtones and combinations of the molecular vibrational modes are found. 
These features are weaker than the fundamental modes (e.g., Gerakines et al. 2005). but 
there is a higher flux of solar radiation in the near-IR. Most work in literature uses the 
wavelengths of absorption features in pure, unprocessed ices in comparison to observed data 
in order to deduce the general composition of icy surfaces (e.g., Roush 2001), however recent 
experiments use ice mixtures (e.g., H 2 O+NH 3 by Moore et al. 2007 and N 2 +CH 4 bj Brunetto 
et al. 2008). There is a need for laboratory mixtures in order to make accurate comparisons 
due to the fact that mathematical addition of individual spectra does not take into account 
the interactions of ice components when present in an intimate mixture. 

Studies of interstellar ices have extended into the near-IR spectral region as well. An 
example includes observations made with SpeX on the Infrared Telescope Facility (IRTF) 
telescope (Vacca et al. 2004). The J ames Webb Space Telescope (JWST) Near Infrared 
Spectrometer (NIRSpec) will perform observations from 0.6 to 5 pm at spectral resolutions 
from 100 to 3000 (Rauscher et al. 2004), and will call for an increase in near-IR data. With 
our data, calculations of ice abundances in star-forming regions could be accomplished. 

A vital aspect of understanding ices on cryogenic surfaces in the ISM is the accurate 
determination of the abundances of the ice components. This relies on accurate laboratory 
measurements of physical parameters called band strengths. For example, the column density 
of a molecule in the solid phase ( N , in units of cm -2 ) as determined from astronomical 
observations of transmitted starlight is commonly calculated from the optical depth spectrum 
t(v) using the formula 

/. r(P) d.u 

N = ^ , ( 1 ) 

where the A l is the band strength of molecular absorption feature i (in units of cm molecule" 1 ), 
and the integration is taken over the wavenumber range of feature i. 

From Equation (1) it is clear that the accuracy of the compositional information that 
may be derived from the astronomical data t{u) is closely tied to the accuracy of the band 
strength A { . Additionally, one cannot, in general, apply the same band strength to all ice 
mixtures containing the same molecular absorption feature, but it must be determined for 
each mixture individually as this parameter is dependent upon the specific physical conditions 
of the local molecular environment (neighboring molecules, crystalline structure, tempera- 
ture, etc.; see, e.g., Gerakines et al. 1995 or Quirico & Schmitt 1997). The determination of 
surface ice composition from a reflectance measurement also depends on the strengths of ice 
absorption features (or more precisely, the optical constants of the ice), as these parameters 
are required as input to surface scattering models used for interpretation. 

The IR absorption features of most molecules of astrophysical interest have band strengths 
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from about 10 -20 to 10“ 16 cm molecule -1 , where the strongest absorption features generally 
lie in the mid-IR. Our laboratory group has already published a study of the near-IR band 
strengths of pure CO, CO 2 , C 3 0 2 , H 2 0, NH 3: CH 4 , and CH 3 OH ice films at low temperature 
(Gerakines et al. 2005). This is a continuation of that work, using the same techniques, to 
the near-IR band strengths of selected mixed ices. Following the methods of Gerakines et al. 
(2005), here we present the near-IR band strengths of CO, C0 2 , CH 4 . and NH 3 as diluted 
in either N 2 or H 2 0, in the mixing ratio of ~5:1 at a temperature of 5K. 


2. Experimental Methods 

The experimental setup and ice creation techniques are identical to those described 
by Gerakines et al. (2005) and similar to those in previous studies by, e.g., Hudgins et al. 
(1993) or Gerakines et al. (1995). To summarize, the sample gases were deposited onto an 
IR- transparent substrate (Csl) at low temperature (~ 5K) inside a high-vacuum sample 
chamber, and the resulting films were analyzed by FTIR spectroscopy. The vacuum system 
pressure was approximately 3 x 10 6 torr at room temperature before cooling the substrate. 
The substrate is attached to the end of a rotatable cold finger that may be cooled by a 
closed-cycle helium refrigerator (Air Products Displex 202N sj'stem). The temperature is 
continuously monitored by a chromel-Au thermocouple and is adjustable by a resistive heater 
element from 5 K up to room temperature. 

The gases to be studied are mixed in the desired mixing ratio in a gas bulb attached 
to a separate manifold then transferred to the high-vacuum system. The sample gases are 
released within the vacuum chamber at a position about 1 cm in front of the substrate. Gas 
how into the system is controlled by a needle valve, and flow rates are adjusted by this valve 
to achieve film growth rates of about 1-5 //m hr -1 . Two of the external ports of the vacuum 
chamber hold IR transparent windows made of KBr, and the entire chamber is situated 
inside the sample compartment of a ThermoMattson Infinity Gold FTIR spectrometer such 
that the IR beam passes at normal incidence through the KBr external ports and the Csl 
substrate. 

The spectrometer optics is capable of automatic switching between the configurations 
of sources, detectors, and beamsplitters without breaking the dry air purge, allowing for the 
measurement of both near-IR spectra over the range of 10.000 to 3500 cm -1 (1.0 to 2.86 /jra) 
and mid-IR spectra from 4,500 to 400cm -1 (2.22 to 25.0 /mi). All spectra in this study 
were taken at a resolution of 4 cm -1 , which is suitable due to previous studies indicating the 
FWHM of these features was larger than 4 cm -1 (Gerakines et al. 2005). Ice samples created 
by vacuum deposition at ~5 K are expected to be amorphous rather than crystalline, and 
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this v/as supported by the absorption profiles in the resulting IR spectra. The gases used 
and their purities are as follows: N 2 (Matheson, 98+%), CO (Matheson, 99.9%), C0 2 
(Matheson, 99.8%), CH 4 (Matheson, 99+%), NH 3 (Matheson, 99+%), and H 2 0 (distilled 
by freeze-thaw cycles under vacuum). 

By assuming a linear baseline, the integrated absorbance of each feature (in units of 
cm -1 ) was measured and the errors in the feature areas were estimated by integrating across 
the same limits when no feature was apparent above the noise level in the infrared spectrum. 
For narrow features (A v < 10cm _1 ), the errors were found to be about 0.006-0.01 cm -1 . 
This varied in direct proportion to a feature’s width. 

For a given molecule, the relative strengths of two absorption features were determined 
from the relationship between the areas of these bands as the ice was grown. While one could 
use a single ice spectrum and determine the relative strengths of any two features by taking 
the ratio of their areas, since the near-IR features are much weaker than those in the mid-IR 
(~10-100 times smaller; see absorbance scales in Fig. 1 and Fig. 2), errors in their measured 
areas are relatively much higher. Systematic errors involved in measuring the areas of the 
near-IR features are significantly reduced by monitoring the relationship between two given 
absorption features over the course of a slow deposit. 

In the experiments presented in the next section, we created mixtures of CO, C0 2 , CH 4 , 
and NH 3 diluted in N 2 and diluted in H 2 0 with approximate compositions of dilutant:sample 
of 5:1. These mixtures were chosen in order to represent some of the most general cases 
in both the ISM or the outer Solar System. Dilution in H 2 0 is meant to represent the 
environment of interstellar icy grain mantles or icy moons, such as Europa or Enceladus 
(Brown 2006, Hansen & McCord 2008). Dilution in N 2 is meant to represent the environment 
of icy moons dominated by N 2 , such as Triton (Buie & Grundy 2000, Stansberry 2004). 
Mixture ratios were not meant as a literal representation of astrophysical bodies, rather, 
mixtures were meant to be dilute enough to simulate a N 2 - or II 2 0-dominated environment. 

For each sample molecule, the selected mixture was prepared in the gas manifold as 
described above. After cooling the substrate to ~ 5 K, the gases were slowly bled into the 
vacuum system and the ice film was allowed to grow continuously over the course of several 
hours to days. The film growth was monitored at regular intervals, taking both near- and 
mid-IR spectra throughout the growth period. At each stage, the integrated areas of the 
mid- and near-IR features were measured. When it was determined that most of the near-IR 
absorption features were of sufficient strength to establish a clear trend, the gas flow was 
stopped and the experiment was ended. For each mixture involving H 2 0, it was necessary 
to replenish the gas bulb during the growth in order to deposit a sufficiently thick sample. 
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3. Results 

The spectra of ices containing solid CO, C0 2 , CH 4 , and NH 3 diluted in N 2 in a ~5:1 
mixture ratio and of ices containing solid CO, C0 2 , CH 4 , and NH 3 diluted in H 2 0 in a ~5:1 
mixture ratio were collected in the near-IR region from 10,000 to 3,500 cm -1 (1 .0-2.9 pm) 
and the mid-IR region from 4,500 to 400 cm -1 (2.2-25 pm) during the slow growth of films 
at ~5 K. Figure 1 displays the spectra of the N 2 -dominated samples and Figure 2 displays 
the spectra of the H 2 0-dominated samples. If possible, the band strengths of any near-IR 
absorption features observed in the spectra were determined using the method described 
below. 

For optically thin ice absorptions, depositing a certain number of molecules should cause 
a linear increase the areas under all absorption features, where the rate of increase in area 
versus column density is equal to that feature’s band strength- see Eq. ( 1 ). Hence, for a 
given molecule, the trend in one particular feature’s integrated area versus that of another 
feature of the same molecule should initially be a straight line whose slope corresponds to the 
ratio of the band strengths of the chosen features. In our case, we have chosen to compare 
the integrated areas of the observed near-IR features of each of the molecules under study 
to the areas of their much stronger mid-IR features. The mid-IR features used for these 
comparisons and their previously determined band strengths are listed in Table 1. Figures 3 
- 6 display the near-IR regions of the spectra where absorptions were observed as well as 
the plots of the near-IR feature areas versus mid-IR feature areas for the N 2 -dominated 
ice mixtures. Figures 7-10 display the corresponding spectra and band area plots for 
the H 2 0-dominated ice mixtures. For peaks where there is an underlining feature of H 2 0 
(see Figure 8 for example), we drew a linear baseline underneath the feature to remove the 
underlying absorption due to II 2 0. 

Once the slope of the linear trend in the near-IR area versus mid-IR area is determined, 
it may be multiplied by the known band strength of the mid-IR feature (from Table 1) to 
obtain the near-IR band strength. Mid-IR band strength data exist for mixtures similar to 
the H 2 0-ice mixtures and the N 2 + CH 4 mixture studied here (Brunetto et al. 2008; Kerkhof 
et al. 1999; Gerakines et al. 1995) and those results were used to determine our near-IR band 
strength values. For GO, C0 2 , and NH 3 in N 2 -dominated mixtures, no such studies have 
been reported in the mid-IR. However, the studies of the molecules in H 2 0 have shown that 
the mid-IR band strengths change by less than 20% relative to their pure-ice values (e.g. 
Gerakines et al. 1995). Since N 2 is a non-polar molecule, it can be expected that its effects 
would be even less pronounced. This is consistent with data for CH 4 (Kerkhof et al. 1999; 
Brunetto et al. 2008). Hence, where data for N 2 -dominated mixtures are unavailable, we have 
used the mid-IR band strengths for pure ices, with the caveat that this may introduce an error 
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of approximately 20 % (but this is within the typical error in our measurements; see results 
below). Observed feature parameters (positions, widths, and vibrational assignments), slopes 
obtained for the linear trends in the near-IR vs. mid-IR areas, and calculated band strengths 
using this method are given in Table 2 for N 2 -dominated mixtures and in Table 3 for H 2 0- 
dominated mixtures. The listed uncertainties in near-IR band strengths were obtained by 
multiplying the standard deviations in the best-fitting slopes bv the mid-IR band strengths 
and do not include the systematic errors associated with the band strength used for scaling. 

For some of the molecules (CO, C0 2 , and CH 4 ), the scaling process was complicated 
by the fact that some mid-IR absorptions became saturated (with peak absorbance values 
greater than 1.0) early in the ice film growth. The curves of growth were highly non-linear. 
Thus, weaker features needed to be used in order to produce linear trends. For the CO ice 
deposits the absorption feature of 13 CO at 2092 cm -1 (4.780 gm) was used to scale the 12 CO 
and 13 CO near-IR features. The strength of 1.5xl0 -19 cm per 12 CO molecule was used 
for the 2092 cm -1 (4.780 gm) feature, taking the band strength of the 13 CO feature from 
Gerakines et al. (1995) and scaling by the terrestrial ratio of 12 C/ 13 C of 87. For N 2 + C0 2 
(5:1), the combination mode at 3711cm -1 (2.695 //m) was used. For the H 2 0 + C0 2 (5:1) 
ice deposit, the 2280 cm -1 feature of 13 C0 2 was used in the same manner. For the H 2 0 + 
CH 4 (5:1) ice the mid-IR absorption feature at 1306 cm -1 (7.657 //m) was used for scaling 
(strengths are listed in Table 1). 

In the N 2 + CH4 ice, the feature due to the CH 4 , the i^ + zq vibration mode at 2825 cm -1 
(3.540 ^m), was used, assuming a band strength of A = 2.8xl0 -19 cm molecule -1 (corre- 
sponding to the 2815 cm -1 feature in pure CH 4 ice as determined by Brunetto et al. 2008). 
Note, in Gerakines et al. 2005, a scaling error caused each of the pure CH 4 results to be off 
by a factor of 7 and was corrected for future work. 

A small level of contamination due to C0 2 (~2 %, as determined from the size of the 
absorption) was also observed within the N 2 ice mixtures, as is seen in Figure 1. Additionally, 
a very small amount of H 2 0 was seen in the N 2 ice mixture (< 1 %) . The value is small enough 
to not be seen in the Figure 1. due to the short time scale required for N 2 -dominated mixtures 
versus the H 2 0-dominated mixtures, which required changing of gas bulbs and additional 
pumping time between scans (as described in the Experimental Methods Section). Near-IR 
features occurring at ~7300 cm -1 and ~5200cm -1 were observed in many H 2 0-dominated 
spectra (e.g., in Figure 2). These features were determined to be caused by small amounts 
of C0 2 and N 2 contamination within the H 2 0 used to create the gas mixtures. This was not 
expected to significantly affect the band strength results. 



4. Comparisons with previous studies 


In this section, we compare these results to those found in the study of pure ices by 
Gerakines et al. (2005), and to other works in the literature. 


4 . 1 . N 2 + CH 4 

For CII4 diluted in N 2 , no visible features were detected for the (i/ 2 + 2v$ + 1/4) and 3 iaj 
vibrational modes near 8780 and 9060 cm -1 , however this may have been due to noise in the 
spectrum. All other vibrational modes showed an increase in peak wavenumber relative to 
pure CH 4 (ranging from 10 cm -1 to 23 cm -1 respectively). All of the features were broader 
than in pure CH 4 (by 2 cm' 1 to 11 cm -1 ), with the exception of the (i/ 4 -I- 1/ 4 ) mode, which 
saw no broadening or sharpening of the peak, and the [i/j] + [10] + [n 4 ] + [1/4] and the (2i/ 3 -f-2i/ 4 ) 
modes, which saw sharpening of the peaks and a decrease in the FWHM by 4 cm -1 and 7 
cm -1 respectively. Band strengths for the (1/3 + U4), (u 2 + u s ), 2 i/ 3 , and {y 2 + 2 i/ 3 ) features (at 
4311, 4539, 6000, and 7503 cm -1 respectively) measured were lower than those in the pure 
ice by a factor of 3 - 13. The (i/ 2 -f ^ 4 ), + i/ 4 ), ( v 3 + 2i/ 4 ), and (i/ 2 + 1/3 + i/ 4 ) features (at 

4125, 4222, 5576, and 5811 cm -1 respectively) saw little to no change in the band strengths 
in comparison to the pure ice (< 10 %). Due to the non-linear behavior of the growth curve 
at short deposition times, the strength of the [t/j] + [z/ 4 ] + [zvj] + [2/4] feature at 8422 cm -1 
could not be determined with accuracy. 

Near-IR spectra of N 2 -I- CH 4 in (1:1) and (10:1) ratios were reported by Brunetto et al. 
(2008), and our band strengths are comparable. For all features, our 5:1 mixtures indicated 
slightly larger peak shifts relative to the pure ice of an additional 3 14 cm -1 . 


4.2. N 2 + C0 2 

For C0 2 diluted in N 2 , all expected near-IR vibrational modes were observed and showed 
a positive in peak position relative to the pure C0 2 ranging from 3 cm -1 to 11 cm -1 , with 
the exception of the (2i/ 4 + 2v 2 + v 2 ) mode and the mode, which showed a decreases in 
position of 1 and 2 cm -1 respectively. A slight broadening relative to pure C0 2 occurred 
in most features, ranging from 1 to 2 cm -1 , with the exception of the (2^ + 2v 2 + v 3 ) and 
(2l>i + j/ 3 ) modes, which showed no broadening or sharpening, and the (4u 3 + u 2 ) feature, 
which showed a sharpening by 6 cm -1 . All measurable band strengths showed an increase 
by 1.3 to 2 times those of pure C0 2 . Only small broadening and shifting occurred due to 
the nature of the apolar ices in the mixture (C0 2 and other apolar ices are not typically self- 



interactive). The increase in band strength may have been due to the slight contamination 
mentioned in Section 3. 


4.3. N 2 + CO 

For CO diluted in N 2 , all expected near-IR. vibrational modes were observed and showed 
little change in peak position relative to pure CO (+1 cm -1 for all modes, except the 2v 13 CO, 
which showed no change in position). Broadening relative to pure CO occurred in each 
feature ranging between 1 and 5 cm -1 . The only measurable band strength was that of the 
2 v 12 CO vibrational mode. The band strength was very comparable to the pure ice feature’s 
band strength (1.1±0.1) x 10 -19 cm molec -1 versus (1.6±0.1) x 10 -19 cm molec -0 for the 
pure ice). As was seen with C0 2: CO saw only minimal broadening and shifts and almost 
no change in the band strength due to the nature of two apolar ices in mixture, indicating 
almost no change in the interaction of the molecules. 


4.4. N 2 -f- NH3 

For NH3 diluted in N 2 . all expected near-IR vibrational modes were observed, with large 
shifts in the peak positions relative to pure NH 3 . Decreases in position occurred in both the 
(yj + y 2 + ^4) and the {y x + y 2 ) modes (of 65 and 30 cm -1 respectively), while increases were 
seen in the peak positions of the 2u x and (u x + 1/4) features (of 83 and 49 cm -1 ). All features 
showed a sharpening of peaks relative to pure NH 3 with decreases of 55 cm -1 to 97 cm -1 . 
There was a decrease in the band strengths for the (u x + y 2 ) and (u x + i/ 4 ) features by a 
factor of 2 to 4 in comparison to those of pure NH 3 and an increase in the strength of the 
2^i mode feature by a factor of 2. When diluted with N 2 , the mixture has fewer NH3 NH 3 
interactions, likely resulting in the sharpening of the peaks. 


4.5. H 2 0 + CH 4 

For CH 4 diluted in H 2 0, no visible absorption features were detected for the (y 2 + 2 y 4 ) 
vibrational mode near 4115 cm -1 . All other vibrational modes showed an increase in peak 
wavenumber relative to pure CH 4 (ranging from 4 to 213 cm -1 respectively), with the 
exception of the (u 2 + v$) and modes at 4528 and 9039 cm -1 (the peak position of 

the 3i/ a mode decreased by 21 cm -1 from the pure ice value). Most of the CH 4 absorption 
features in the H 2 0 mixtures were broader than their corresponding pure-ice rallies, with 
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increases in FWHM from 1 to 19 cm -1 , with the exception of the [z/j] + [u\] + [z/ 4 ] 4- [i/ 4 ], 
(2 z/ 3 + 2^ 4 ), and 3z/ 3 modes (at 8569, 8612, and 9039 cm -1 ), which all showed a sharpening 
of the features with respect to their pure-ice values by 6 — 35 cm -1 . Due to the non-linear 
behavior of the growth curves even at short deposition times, the strengths of only two near- 
IR features could be determined with accuracy, as listed in Table 3, and each were found to 
be weaker than in the pure ice. 

A H 2 0 -I- CH 4 mixture was studied at 15 K in two dilution ratios (20:1 and 3:1) by 
Bernstein et al. (2006) from 1.5 to 25 /*m (6700 to 400 cm -1 ). Generally, our measured 
widths and band strengths agree with those reported for the (1^ + z/ 4 ), (z/ 3 4- v 4 ), (u 2 + z/ 3 ), 
(v 2 + v 3 4- Vi), and 2^3 modes. For all features, our 5:1 mixtures displayed slightly higher 
peak shifts than Bernstein et al. (2006) relative to the pure ice of an additional 4—10 cm -1 . 

Galvez et al. (2009) also studied several mixtures of H 2 0 + CH 4 at 14 K. We find 
that our results agree in widths and band strengths for their H 2 0 4- CH 4 (10:1) and (4:1) ice 
mixtures for the ( v : + v 4 ) and (z/ 3 4- V\) modes. Peak positions were also similar, with only 
a differences of 5 - 8 cm -1 . 


4.6. H 2 0 C0 2 

For C0 2 diluted in H 2 0, no visible features were detected for the {A.v 2 + u 3 ) and (4z/ 3 ^z/ 2 ) 
vibrational modes near 4840 and 9955 cm -1 . The (i^ + 2z^ 2 + r' 3 ) at 4965 cm -1 , decreased its 
peak positions bv 6 cm -1 respectively in comparison to pure C0 2 value, while the (2i/\ 4- v 3 ) 
mode at 5087 cm -1 saw a slight decrease in position by 1 cm -1 . Larger shifts were seen for 
the (z/i + 4u 2 + vs) mode at 6051 cm -1 (decrease of 163 cm -1 ), the (2z/ 4 + 2v 2 4- i/ 3 ) mode at 
6325 cm -1 (decrease of 16 cm -1 ), and the 3^ 3 mode at 6864 cm -1 (decrease of 108 cm -1 ). 
While broadening occurred in every peak, the most prominent peaks saw the least amount 
of increase in FWHM. The (z/j + 2v 2 + z/ 3 ) and (2v\ + v 3 ) modes at 4965 and 5086 cm -1 saw 
broadening of only 2 and 3 cm -1 respectively, while the (z/j + 4z/ 2 -(- z/ 3 ) , (2z/ 4 + 2z^ 2 + z/ 3 ) . and 
3^3 modes at 6051, 6325, and 6864 cm -1 saw much larger increases in FWHM (93, 71. and 
171 cm -1 respectively). The strengths of only two absorption features could be determined 
with confidence: the (1/4 + 2v 2 + u 3 ) and [2u\ + 2v 2 + u 3 ) modes, which showed a strength of 
2-4 times those of the pure C0 2 band strengths. 

H 2 0 + C0 2 was also studied by Bernstein et al. (2005) at 15 K in two mixtures, 25:1 
and 5:1, from 2 to 25 /im (5000 — 400 cm -1 ). Several near-IR bands were detected, including 
the (ui + 2z / 2 + z/ 3 ) feature, which in the 5: 1 mixture showed a lower band strength 3.4 x 10 -20 
cm molecule -1 compared to our value of (2.3±0. 1) x 10 -19 cm molecule -1 and a larger FWHM 
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(15 cm -1 compared to 9 cm -1 ) than in our results. In the 25:1 mixture of Bernstein et al. 
(2005), the band strength and FWHM of this feature was consistent with those in our 5:1 
mixture. In each case, we measured slightly different peak positions. The (2 i/ 4 + Va) mode 
v'as not detected by Bernstein et al. (2005) in their 25:1 mixture; however, the properties of 
this feature in the 5:1 mixture were similar to those measured here. Bernstein et al. (2005) 
also observed two additional peaks, the 2 and 4 i/ 2 + v$ modes, that were not visible in our 
spectra. This is likely due to overlapping H 2 0 features (see Figure 2 or Figure 8). Another 
reason for this may be that the 2 mode is a forbidden overtone of C0 2 and is extremely 
sensitive to conditions such as concentration and temperature. 


4.7. H 2 0 + CO 

For CO diluted in H 2 0, only the 2u 12 CO and Zv 12 CO features were visible at 4260 
and 6339 cm -1 respectively. Compared to a pure ice, each showed an increase in peak 
position (bv 8 cm -1 and 1 cm” 1 , respectively) and each broadened, by 30 cm -1 and 1 cm -1 
respectively. Only the strength of the 2u 12 CO feature was able to be determined from our 
data ((1.4T0.1) x 10 -19 cm molecule” 1 ), and it was found to be comparable to the strength 
in the pure ice ((1.6 ±0.1) x 10“ 19 cm molecule -1 ). 

4.8. H 2 0 + NH 3 

For the NH 3 diluted in H 2 0, all observed features were seen to shift in peak position 
in comparison to pure NH 3 , with increases ranging from 22 — 82 cm -1 , and to decrease in 
FWHM by 55 — 97 cm -1 . The strengths of all detectable features also decreased from their 
pure ice values by approximately 1.5 times. The strength of the (ui + 2u 2 + mode near 
6136 cm -1 could not be determined from our data. 

H 2 0 + NH 3 (~3:1) at 9K was studied by Moore et al. (2007). Their results included 
peak positions of the (ui ± u 2 ) and (u~ ± ± i/ 4 ) vibrational modes. For the (ui ± i/ 2 ) feature, 
Moore et al. (2007) reported a peak position of 4478 cm -1 , which is further (by 7 cm -1 ) 
from the pure ice value than our results indicate. The ( v i + i/ 4 ) mode also had a larger 
change from the pure ice value than seen in our results (by 2 cm” 1 ). Moore et al. (2007) 
did not calculate bands strengths. 
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5. Discussion and Implications 

These experimental data can be used for more powerful interpretations of observations 
of astrophysical ices. In particular, the band strengths measured here may be used to more 
accurately determine column densities of ice components in the ISM. As shown in Gerakines 
et al. (2005), using our laboratory data, we can predict the optical depths of near-IR ice 
absorptions for several astrophysical objects. We have estimated the optical depths for sev- 
eral lines of sight in the ISM using the widths and strengths (as well as previously published 
column densities) for some of the near-IR features of our H 2 0-dominated mixtures(e.g., Ger- 
akines et al. 1999; Gibb et al. 2004). These estimates are listed in Table 4 and, due to their 
direct dependence on band strengths, show weakening in comparison to the results predicted 
in Gerakines et al. (2005). The estimates here are more appropriate to interstellar ices than 
our pure ices analysis in Gerakines et al. (2005), since interstellar ices are dominated by 
water ice. However they may not fully reflect the most realistic cases for these objects, as 
some are located in highly energetic environments, such as NGC 7538 IRS 9, GL 7009S, 
and GL 2136. which are all in the vicinity of young stellar objects (YSOs) and therefore 
undergoing various degrees of heating and/or energetic processing (Gibb et al. 2000). 

N 2 -dominated ice spectra may be used in comparison with observed planetary ices in- 
stead of pure ices to deduce information about the general composition of the surface ice. 
An example is Fig. 2 in Quirico et al. (1999) in which the spectrum of Triton is compared 
with laboratory transmission spectrum of pure ices CH 4 , CO, N 2 , C0 2 , and H 2 0. Similar 
to Triton, Pluto has been recently investigated by Buie et al. (2010) and has shown signs 
of N 2 ice, however, the behavior (i.e. seasonal effects) of ices on both planetary surfaces 
remains poorly understood. With the upcoming New Horizons mission in 2015, the data 
presented here may be useful for interpretation and further understanding of the general 
composition of the planetary surface. However, the chemical compositions of icy planetary 
bodies are affected by particle radiation or high-energy photons, resulting in the synthesis of 
new molecular species which can then be ejected into the gas phase by various mechanisms 
(Bringa & Johnson 2004). Studying the effect of UV radiolysis on these ices in addition to 
placing them in molecular mixture would lead to more realistic choices for comparison to 
observational data. 


6. Summary and Future Work 

In this paper, we have collected spectra of N 2 mixed with CO, C0 2 , CH 4 , and NH 3 
in 5:1 ratios and of H 2 0 mixed with CO, C0 2 , CH 4 . and NH 3 in 5:1 ratios in the Near- 
IR region from 10,000-3,500 cm _1 (1.0-2.9/rm) and the mid-IR region from 4500-400 cm -1 
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(2.2-25 /xrn) during the slow growth of films at ~5 K and have used these to the determine 
near-IR band strengths listed in Table 2 and in Table 3. Future work will analyze the effect 
of UV radiolysis on both the pure ices previously discussed in Gerakines et al. (2005) and 
the ice mixtures described in this paper. We will also investigate the effects, if any, of ice 
composition on the near-IR band strengths as well as the differences between the near-IR 
band strengths of crystalline and amorphous ices. 
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Table 1. Mid-Infrared Features and Band Strengths Used in Strength Determinations 


Molecule 

v [cm x ] 

A [/xm] 

A [cm molecule : ] 

Reference 

Pure CILj 

1306 

7.657 

7.0 x 10- 18 

1 


2815 

3.552 

2.8 x 10- 19 

2 

CH 4 in H 2 0+ CH 4 (1.03:1) 

1302 

7.680 

8.1 x 10" 18 

3 

N 2 + CH 4 (10:1) 

2820 

3.546 

3.4 x 10 -19 

2 

Pure C0 2 a 

2283 

6.639 

7.8 x 10- 17 

4 


3708 

2.697 

1.4 x 10~ 18 

4 

C0 2 in H 2 0+ C0 2 (1.6:1) 

2279 

6.639 

7.6 x 10- 17 

4 

Pure CO a 

2092 

4.780 

1.5 x 10 19 b 

5 

CO in H 2 0+ CO (2.1:1) 

2092 

4.780 

1.2 x 10- 19 b 

4 

Pure NH 3 a 

1070 

9.346 

1.2 x 10- 17 

5 

NH 3 in H 2 0+ NH 3 (11:1) 

1070 

9.346 

1.3 x 10- 17 

3 


a Used for N 2 mixture. See text for farther discussion. 

b Although this is a feature of 13 CO, its band strength is expressed in units of cm per 12 CO 
molecule. 


Note. — (1) Boogert et al. (1997); (2) Brunetto et al. (2008); (3) Kerkhof et al. (1999); 
(4) Gerakines et al. (1995); (5) d’Henaecourt & Allamandola (1986). 
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Table 2. Observed Near-IR Features and Band Strengths of N 2 -dominated Ice Mixtures at 

5K 


Peak Position 
v [cm -1 ] A [/jm] 

FWHM, tAv 
[cm -1 ] 

Vibration 

Mode 

Relative Growth a 
Rate 

Strength, A 
[cm molecule -1 ] 

N 2 +CH 4 (5:1) 

4125 2.425 

14 

1/2 4- 21/4 

(3.9 ±0.3) x 10 -2 

(1.3 ±0.1) x 10 -20 

4222 

2.374 

11 

Vl + 2/ 4 

(79 ± 4) x 10 -2 

(2.7 ±0.1) x 10 -19 

4311 

2.320 

18 

2/3 + 2 ^ 

(230 ± 4) x 10 -2 

(7.9 ±0.1) x 10 -19 

4539 

2.203 

15 

1*2 ± 2/3 

(27 ±0.4) x 10 -2 

(9.2 ±0.1) x 10 -2 ° 

5576 

1.793 

15 

v 3 4- 2j / 4 

(3.6 ± 0.4) x 10 -2 

( 1.2 ± 0 . 1 ) x 10 -2 ° 

5811 

1.721 

15 

^2 t V 3 + ^4 

(7.1 ± 0.4) x 10 -2 

(2.4 ±0.1) x 10 -2 ° 

6000 

1.667 

18 

2 l /3 

(11 ±0.4) x 10 -2 

(3.6 ±0.1) x 10 -2 ° 

7503 

1.333 

18 

14 2 + 2 l /3 

(20 ± 0.4) x 10 -3 

( 6.8 ± 0 . 1 ) x 10 -21 

8422 

1.187 

26 

M + [td] + M + M 

... 


8611 

1.161 

35 

2^3 + 2 1/ 4 

(5.4 ± 0.4) x 10 -2 

( 1.8 ± 0 . 1 ) x 10 -2 ° 

N 2 +CO 2 (5:1) 
4840 2.066 

7 

4i / 2 + 2/3 

(9.3 ± 0.7) x 10 -3 

(1.3 ±0.1) x 10 -2 ° 

4975 

2.010 

7 

V\ + 2^2 + 2/3 

(59 ±0.7) x 10 -3 

( 8.2 ± 0 . 1 ) x 10 -2 ° 

5090 

1.965 

7 

2 i^i + 2/3 

(25 ± 0.7) x 10 -3 

(3.5 ±0.1) x 10 -2 ° 

6218 

1.608 

8 

V\ + 4z/2 + ^3 



6340 

1.577 

9 

2 V\ + 2 i *2 + V 3 



6970 

1.435 

10 

3^3 



9955 

1.005 

6 

4 i / 3 + 2/ 2 (?) 

... 


N 2 +CO (5:1) 

4159 2.404 

7 

2 u i3 CO 



4253 

2.351 

6 

2 i / 12 CO 

(73 ± 7) x 10 -2 

(1.1 ±0.1) x 10 -19 

6339 

1.578 

9 

3i/ 12 CO 

... 


8505 

1.176 

9 

? 

... 


N 2 +NII 3 (5:1) 
4444 2.250 

22 

Vi 4 - V 2 (?) 

(18 ± 0 . 8 ) x 10 -3 

(2.2 ±0.1) x 10 -19 

5042 

1.983 

13 

4 - 2/4 (?) 

(33 ± 0.8) x 10 -3 

(4.0 ±0.1) x 10 -19 

6034 

1.657 

- 

4- V-2 4 - 2/4 (?) 



6598 

1.516 

21 

22/1 (?) 

(7.2 ± 0.2) x 10 -3 

(8.6 ± 0 .2) x 10 -2 ° 


a The slope of the linear trend in the near-IR vs. mid-IR band areas shown in Figures 3-6. 


Note. — Vibrational mode assignments, except those in NH 3 mixtures, from Calvani et al. 1992 and 
Quirico & Schmitt 1997; NH 3 assignments are uncertain, however, are similar to those made in Moore et al. 
(2007). 
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Table 3. Observed Near-IR Features and Band Strengths in H 2 0-dominated Ice Mixtures 

at 5 K 


Peak Position 
v [cm -1 ] A [/nnj 

FWHM, \v 
[cm -1 ] 

Vibration 

Mode 

Relative Growth a 
Rate 

Strength, A 
[cm molecule -1 ] 

H 2 O+CH 4 (5:1) 
4209 2.376 

17 

V' A +V 4 

(3.6 ±0.1) x 10 -2 

(2.9 ±0.1) x 10 - 19 

4308 2.321 

16 

V 3 + v 4 

(7.0 ± 0.3) x 10 “ 2 

(5.7 ± 0.2) x 10 -19 

4528 2.208 

10 

V 3 + V 3 

... 


5640 1.773 

8 

1/3 + 2 l / 4 

... 

... 

5804 1.723 

23 

V 2 + V 3 + V 4 


... 

5995 1.668 

11 

2 i / 3 


... 

7700 1.299 

26 

Vi + 2 1/3 



8569 1.167 

22 

[^i] + [th] + [^ 4 ] + [^ 4 ] 

... 


8612 1.161 

7 

2 i / 3 + 2 1/4 

... 


8962 1.116 

36 

V 2 + 2 1/ 3 + 1/4 



9039 1.106 

22 

3i / 3 



H 2 O+CO 2 (5:1) 

4965 2.014 

9 

1^1 + 2 l >2 + V 3 

(3.1 ±0.1) x 10 -3 

(2.3 ±0.1) x 10 -19 

5086 1.966 

10 

2 i/i + v 3 

(10 + 0.5) x 10 -4 

(7.6 ± 0.4) x 10 -2 ° 

6051 1.653 

100 

vi + 4i / 2 + V 3 

... 

... 

6325 1.581 

78 

2 v\ + 2 t / 2 + IS 3 


... 

6864 1.457 

179 

3^3 

... 


H 2 O+CO (5:1) 

4260 2.347 

33 

2i/ 12 CO 

1.2 ±0.7 

(1.4 ±0.8) x 10 -19 

6339 1.575 

5 

3i / 12 CO 



H 2 0+NH 3 (5:1) 

4518 2.213 

22 

"1 + ^2 (?) 

(44 ± 0.8) x 10 -3 

(5.7 ±0.1) x 10 “ 19 

5015 1.994 

13 

v 1 + Vi (?) 

(38 + 0.8) x 10 -3 

(4.9 ±0.1) x 10 -19 

6136 1.630 

* 

Vi +v 2 + v 4 (?) 

... 

... 

6597 1.516 

21 

2^1 (?) 

(27 ± 0.8) x 10 -4 

(3.9 ±0.1) x IQ -20 


a The slope of the linear trend in the near-IR vs. mid-IR band areas shown in Figures 7-10. 


Note. — Vibrational mode assignments, except those in NH 3 mixtures, from Calvani et al. 1992 and 
Quirico & Schmitt 1997; NH 3 assignments are uncertain, however, are similar to those made in Moore et al. 
(2007). 
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Table 4. 


Column Densities* 1 


and Predicted Near-IR Optical Depths 1 * of Interstellar Ices 
Diluted in H 2 O 



W33A 

NGC 7538 IRS9 

Elias 16 

Sgr A* 

GL 7009S 

GL 2136 

ch 4 

N 

1.7 x 10 17 

1.5 x 10 17 


2.6xl0 16 

4.8xl0 17 


r( 2.376 H 

0.003 

0.002 

— 

< 0.001 

0.008 

— 

r(2.321/nn) 

0.006 

0.006 


< 0.001 

0.01 

— 

C0 2 

N 

1.4 x 10 18 

1.7 x 10 18 

4.5 x 10 17 

1.8 x 10 17 

2.5 x 10 18 

8.0 x 10 17 

r(2.014/xm) 

0.04 

0.05 

0.01 

0.005 

0.07 

0.02 

r(1.966 fim) 

0.01 

0.01 

0.004 

.001 

0.02 

0.006 

CO 

N 

8.8 x 10 17 

1.2 x 10 18 

6.3 x 10 17 

< 2 x 10 17 

1.8 x 10 18 

1.0 x 10 17 

r( 2.347 /an) 

0.004 

0.005 

0.002 

<0.001 

0.008 

<0.001 

nh 3 

N 

1.7 x 10 18 

9.8 x 10 17 

< 2.3 x 10 17 

2.6 - 3.9 x 10 17 


— 

r( 2.213 /xm) 

0.04 

0.02 

< 0.006 

0.006 - 0.01 

— 


r(1.994 fi m) 

0.06 

0.03 

< 0.009 

0.01 - 0.02 


— 


a Column densities ( N ) have units of cm 2 ; Gibb et al. (2004). 

b Optical depths have been calculated from the listed column densities, using width and strength data 
listed in Table 3. 
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Fig. 1. — Combined near- and mid-IR spectra (10,000-400 cm -1 , 1.0-25 /xm) of ice films at 
T ~ 5 K for the N 2 -dominated mixtures. The sinusoidal baseline behaviors (strongest in the 
N 2 + C0 2 and N 2 + CO samples) are due to interference of the IR beam within the samples. 
Some features due to minor contamination by H 2 0 and C0 2 are also visible (see text for a 
full discussion). 

Fig. 2.— Combined near- and mid-IR spectra (10.000-400 cm -1 , 1.0-25 /xm) of ice films at 
T ~ 5K for the H 2 0-dominated mixtures. The sinusoidal baseline behaviors are due to 
interference of the IR beam within the samples. Some features due to minor contamination 
by C0 2 are also visible (see text for a full discussion). 

Fig. 3. — Top: Selected near-IR absorbance spectra of a N 2 + CH 4 (5:1) ice mixture at 
T ~ 5K for the most prominent features studied. Bottom: Integrated areas (in cm -1 ) of 
CH 4 features vs. the area (in cm -1 ) of the 2825 cm -1 (3.540 /x m) feature during deposition 
at T ~ 5 K. The solid lines are linear fits to the solid symbols, and data points shown as 
empty symbols were omitted from the fits (see explanation in text). 

Fig. 4. — Top: Selected near-IR absorbance spectrum of a N 2 -f C0 2 (5:1) ice mixture at 
~ 5K for the most prominent features studied. Bottom: Integrated areas (in cm -1 ) of C0 2 
features vs. the area (in cm -1 ) of the 3711 cm -1 (2.695 /xm) feature during deposition at 
r~> 5 K. The solid lines are linear fits to the solid symbols, and data points shown as empty 
symbols were omitted from the fits (see explanation in text). 

Fig. 5. — Top: Selected near-IR absorbance spectrum of a N 2 + CO (5:1) ice mixture at 
~ 5 K for the most prominent features studied. Bottom: Integrated areas (in cm" 1 ) of CO 
features vs. the area (in cm" 1 ) of the 2092 cm -1 (4.780 /xm) feature during deposition at 
~ 5 K. The solid lines are linear fits to the solid symbols. 

Fig. 6 . — Top: Selected near-IR absorbance spectra of a N 2 + NH 3 (5:1) ice mixture at 
~ 5K for the most prominent features studied. Bottom: Integrated areas (in cm -1 ) of NH 3 
features vs. the area (in cm -1 ) of the 1146 cm " 1 (8.726 /xm) feature during deposition at 
~ 5 K. The solid lines are linear fit to the solid symbols, and data points shown as empty 
symbols were omitted from the fits (see explanation in text). 

Fig. 7. — Top: Selected near-IR absorbance spectrum of a H 2 0 -t- CII 4 (5:1) ice mixture at 
-5K for the prominent most features studied. Bottom: Integrated areas (in cm -1 ) of CH 4 
features vs. the area (in cm -1 ) of the 1306 cm -1 (7.657 /xm) feature during deposition at 
— 5 K. The solid lines are linear fits to the solid symbols. 
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Fig. 8. — Top: Selected near-IR absorbance spectrum of a H 2 0 + C0 2 (5:1) ice mixture at 
~ 5K for the most prominent features studied. Bottom: Integrated areas (in cm' 1 ) of C0 2 
features vs. the area (in cm -1 ) of the 2280 cm' 1 (4.386 jum) feature during deposition at 
5 K. The solid lines are linear fits to the solid symbols, and data points shown as empty 
symbols were omitted from the fits (see explanation in text). 

Fig. 9. — Top: Selected near-IR absorbance spectrum of a H 2 0 + CO (5:1) ice mixture at 
~ 5K for the most prominent features studied. Bottom: Integrated areas (in cm -1 ) of CO 
features vs. the area (in cm -1 ) of the 2092 cm -1 (4.386 pm) feature during deposition at 
~ 5 K. The solid lines are linear fits to the solid symbols. 

Fig. 10. — Top: Selected near-IR absorbance spectra of a H 2 0 + NH 3 (5:1) ice mixture at 
~ 5 K for the most prominent features studied. Bottom: Integrated areas (in cm -1 ) of NH 3 
features vs. the area (in cm -1 ) of the 1110 cm' 1 (9.009 pm) feature during deposition at 
~ 5 K. The solid lines are linear fits to the solid symbols, and data points shown as empty 
symbols were omitted from the fits (see explanation in text). 
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FIGURE 4 
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FIGURE 5 
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FIGURE 7 
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FIGURE 8 
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FIGURE 9 
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